Nonequilibrium Critical Dynamics of the 2D XY model 
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The nonequilibrium critical dynamics of the 2D XY model is investigated numerically through Monte 
Carlo simulations and analytically in the spin-wave approximation. We focus in particular on the 
behaviour of the two-time response and correlation functions and show that the ageing dynamics 
depends on the initial conditions. The presence of critical fluctuations leads to non-trivial violations 
of the fluctuation-dissipation theorem apparently reminiscent of the three dimensional Edwards- 
Anderson spin glass model. We compute for this reason the finite-size overlap probability distribution 
function and find that it is related to the finite-time fluctuation-dissipation ratio obtained in the 
out of equilibrium dynamics, provided that the temperature is not very low. 
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I. INTRODUCTION 

One of the most interesting problems in the study of 
glassy systems [jjj is the characterization of the nonequi- 
librium dynamics in terms of equilibrium phase space 
structure. In nonequilibrium dynamics temporal correla- 
tions do not show time translational invariance. In the 
so-called ageing regime the correlation function C(t,t w ) 
of any observable is a two-time function depending on 
both a waiting time t w and t, the total time. C(t,t w ) 
and its conjugate linear response R(t,t w ) are related in 
a suitable long-time limit by a generalized form of the 
fluctuation-dissipation-theorem (FDT) M~m, 



R(t, t w ) 



x(C(t,t w )) dC(t,t w ) 



k B T 



(1) 



The factor x(C(t,t w )), the fluctuation-dissipation ratio 
(FDR), is equal to 1 at equilibrium and its departure 
from this value characterizes the ageing dynamics of a 
system out of equilibrium. A simple kind of ageing, 
with x = 0, occurs for example in an Ising ferromag- 
net quenched from a random configuration into the or- 
dered phase below the critical temperature ||. Gener- 
ally, glassy systems are characterized by more complex 
behaviour, with the FDR a non-trivial function of t and 

In Ref. km, Franz et al. have shown that in a certain 
class of systems the FDR is closely related to the Parisi 
function which measures the distribution of overlaps be- 
tween pairs of pure equilibrium states . The FDR is 
linked to the Parisi function P(q) through the relation |q] 



x(C) 



dqP(q). 



(2) 



The relevance of such a kind of 'ergodic' result is ob- 
vious as it allows one to predict the long time value 



of physical observables in a non-stationary dynamical 
regime, in terms of static averages computed at equi- 
librium, i.e. using the Gibbs measure. From this equa- 
tion, a complex form for x(C) is related, through P(q), 
to a system with many pure equilibrium states. The re- 
sult x = for ferromagnets quenched below the criti- 
cal point H is also consistent with it, as in this case 
P(q) = S(q — m 2 ). Further, it has been tested numer- 
ically with success in finite-dimensional spin glasses [0. 
However, many interesting situations also arise when the 
technical hypothesis under which the result of Ref. H 
holds are not strictly met. This occurs for instance when 
the asymptotic value of the free-energy is different from 
that at equilibrium p|j9|JTc[. In these examples, viola- 
tions of the fluctuation-dissipation theorem reminiscent 
of glassy systems may indeed appear even in systems with 
trivial Gibbs measures, at least on finite but large time- 
scales jn]]. In this case, one may wonder whether an 
appropriate generalization of (0) can be envisaged. 

It has been recently suggested that non-trivial FDT vi- 
olations appear even in non-disordered systems such as a 
ferromagnet, provided they are at criticality Jl^,|l^]. This 
proposition would be particularly interesting if it led to 
simple models with behaviour analogous to that of glassy 
systems, but for which we have a detailed knowledge of 
the microscopic excitations. In this paper we investigate 
this proposition by studying the nonequilibrium critical 
dynamics of the 2D XY model and its possible relation 
with the static properties. 

At the critical point T = T c , the dynamics 
of a ferromagnet undergoes so-called 'critical slowing 
down' Jl4|,[l5|]. This is a consequence of the fact that at 
T c the equilibrium correlation length £ e q (T e ), and hence 
the largest relaxation time, diverges H,n5[. In par- 
ticular, when a ferromagnet is quenched from a highly 
disordered state to the critical point T = T c |l6|, the 
correlation length increases with the time t elapsed af- 



f 



ter the quench as £(t) ~ t 1 ^, which defines the dy- 
namical exponent z. At time t, only the critical fluc- 
tuations with wavelength smaller than £(t) have equili- 
brated. The magnetization over a length scale is 
given by m^{t) ~ g-P/ v ~ t-Pl vz anc j vanishes at long 
time. In the thermodynamic limit, the system does not 
equilibrate since the equilibrium correlation length is in- 
finite at the critical temperature, £ eq (!T c ) = oo. This is 
rather different from the situation in which a ferromagnet 
is quenched below the critical temperature to T < T c [p.7| . 
Then, well-defined domains with finite and opposite mag- 
netization ±m eq (T) appear and coarsen. At long time t, 
the system enters a scaling regime where the typical size 
L(t) of a domain grows as L(t) ~ i 1 / 2 , where z is the 
dynamic critical exponent defined above fl7j| . This scal- 
ing holds in the regime where £, e q(T) <C L(t) <§C L. Here, 
although equilibrium is never reached in the thermody- 
namic limit, L — > oo, the loss of equilibrium is restricted 
to the domain walls whose density becomes zero at long 
time. The thermodynamics is given by the bulk, which 
is in equilibrium ||. 

The finite-size equilibrium fluctuations of the magne- 
tization to of a critical system extend over all the range 
of length scales between the lattice spacing a and the 
system size L. The probability distribution function 
of to, denoted by Q(to, L), scales then as Q(m 1 L) ~ 
L l3 / u Q{mL f3 / v ), where (3 and v are the usual critical ex- 
ponents associated with the magnetization and the equi- 
librium correlation length, respectively fjs| ]. At the crit- 
ical point the scaling function Q is non-Gaussian and 
markedly asymmetric, with large fluctuations below the 
mean |18|. Hence, for finite-size critical systems a non- 
trivial Parisi function, which for pure ferromagnets is di- 
rectly related to the magnetization distribution, is ex- 
pected to occur. 

In relation (g) , it is understood that the left hand side 
is computed in the limit of large waiting times, after the 
thermodynamic limit is taken. The right hand side is 
computed in an infinite system at equilibrium. In this 
paper, we make the conjecture that FDT violations on a 
finite-time scale t w are governed by the Parisi function 
P(q, L) of a system of finite-size L such that L = £(t w ). 
This amounts to a finite-time, finite-size generalization 
of the relation (0) , namely 



x(C(t,t w )) = 



C(t,t„) 



dqP(q,£(t w )). 



(3) 



The right hand side is now computed by equilibrating a 
finite system. 

To test this proposition, we focus on the static and 
dynamic properties of the 2D XY model below the 
Kosterlitz-Thouless transition, where the system is char- 
acterized by a line of critical points [[l9| . In this respect 
the 2D XY model is particularly interesting, since the 
static distribution of the order parameter fluctuations 
and the non-equilibrium dynamics can be analytically ob- 
tained in the spin-wave approximation. It has also the 



advantage that it is critical over all the low-temperature 
phase and therefore does not require a temperature fine- 
tuning in numerical simulations. Moreover, a similar- 
ity between the 2D XY model and the 3D Edwards- 
Anderson model has often been noticed [^0). Actually, 
the apparent analogy might not be completely acciden- 
tal: it has been proposed that extremal statistics are 
relevant in disordered systems and turbulence |2jJ , while 
Q(to, L) in the 2D XY model is closely related to a gen- 
eralized Gumbel distribution |^2|,^3| which also describes 
the power fluctuations in a confined-turbulence experi- 
ment Q. 

The paper is organized as follows. In the next section, 
we study analytically within the spin-wave approxima- 
tion the behaviour of two-time correlation and response 
functions as well as the FDT violations starting from a 
completely ordered state. In Section III, the same quanti- 
ties are investigated by Monte Carlo simulations for both 
cases of random and ordered initial conditions. This al- 
lows one to go beyond the spin-wave approximation and 
see how the vortex contribution changes the dynamic 
scaling of two-time functions. In Section IV, we compute 
the static fluctuations of the spin glass order parameter 
in a finite size system and test relation (||). Our results 
are finally discussed in Section V. 



II. AGEING DYNAMICS: ANALYTICAL 
RESULTS 

The coarsening dynamics of the 2D XY model has 
been lar gely studied both numerically and theoreti- 
cally [p5|-p3[. These studies have mainly focused on the 
coarsening process itself, rather than on the ageing prop- 
erties (two-time functions and FDT violations). The only 
paper dealing with FDT violations in the XY model is 
Ref. p3| , but the spin-spin autocorrelation function, to- 
gether with its conjugated response function have not 
been studied. We shall fill the gap in this Section, mak- 
ing use of the formalism and results of Refs. [j26|]33| . 



A. Dynamics in the spin-wave approximation 

The relaxational dynamics of the XY model can be 
described by a Langevin dynamics p7?j 



d<p(x, t) 



5(f>(x, t) 



COM)- 



(4) 



In this expression, <j>(x, t) is a two-component order pa- 
rameter, and the free-energy functional F[(j>] is of the 
Ginzburg-Landau type. In the spin-wave approximation 
the free energy functional reads P6J 



P(T) 



(5) 



2 



where p(T) is the spin- wave stiffness and the angular vari- 
able 9 is related to the spin variable by <fi = exp(i6). In 
this approximation the dynamics is described by the fol- 
lowing Langevin equation 



dO(x,t) 
dt 



SF[0] 
~69{x,t) 



+ ((x,t). 



(6) 



In this expression, £ represents the thermal noise. It is 
a random Gaussian variable with mean (C(x, t)) = and 
variance (((x, t)((x', t')) = Anrj(T)p(T)5(x-x')S(t-t'), 
where rj(T) is the usual critical exponent at tempera- 
ture T, linked to the spin-wave stiffness by the relation 
2tt V (T)p(T)=T g|. 

To be complete, initial conditions have to be associ- 
ated with this first order differential equation. Since the 
Hamiltonian (^J) only describes the spin-wave excitations 
of the system, the dynamics (0) cannot be a good rep- 
resentation of the coarsening process following a quench 
from a completely disordered state. In that case vortices 
which are not described by the spin- wave approximation 
are present in the system. We therefore take as initial 
conditions a completely ordered configuration, 



(x,0) = 9 , Vsu. 



(7) 



Similar initial conditions have been studied in Ref. [p3[ , 
but one could generally study quenches from any tem- 
perature below the transition, as was done in Ref. |2(| . 

The dynamics is easily solved by Fourier transforming 
the dynamical equation. This gives 

0(M)= [ dt'e-^P^-^a^t'), (8) 



where 9{k,t) = J d 2 x9(x, t)e~ lk x . The initial condition 
(0), which in Fourier space reads 9(k, 0) = 0, has been 
used. 

The spin-spin autocorrelation function is defined by 

C(t,t w ) = ^ J d 2 x(<j>(x,t) -<t>{x,t w )). (9) 

It can be expressed in terms of the angular variables as 

C(t,t w ) = i J d 2 x{cos[d(x,t)-9(x,t w )}). (10) 

This quantity is easily computed by using the fact that 
6(k,t) is a Gaussian variable, Eq. (g), by developing the 
trigonometric functions and using Wick's theorem 
The correlation reads finally 



C(t,t w ) 



(« 2 + t + OT (T)/4 



(a 2 + 2t)(a 2 + 2t w ) (a 2 + t - t w ) 2 



(11) 



An ultra-violet cutoff a, simulating the lattice spacing has 
been introduced to regularize divergent expressions. This 



is done, for convenience, by adding a factor exp(— k 2 a 2 ) 
in each integration over the Fourier space. 
The conjugate response function is given by 



R(t, t v 



1 

V 



d z x 



6(<Kx,t)) 



Sh(x, tw) 



(12) 



h=0 



where the field h(x,t) is conjugated with <fi(x, t). It may 
also be computed in terms of the angular variables, and 
in a similar way one finds 



R(t, t w ) 



a 4 {a 2 +t + t w ) 2 
(a 2 + 2t)(a 2 + 2t w ) 



n(T)/i 



(13) 



We compute a two-time fluctuation-dissipation ratio 
X(t, t w ) defined as 



X(t, t w ) 



TR(t, t w ) 

dC(t,t m ) ~ 



(14) 



It is given by 
X[t, t w ) — 



(a 2 + t-t w )(t-t w ) 
(a 2 +t + t w )(a 2 + 2t w ) 



(15) 



B. Scaling behaviour 

The dynamic functions C(t,t w ), R(t,t w ) and X(t,t w ) 
may be rewritten in the following scaling forms, for times 
t - t w > a 2 : 



C(t, t w ) = 
R{t,t w ) — 
X{t, t w ) 



1 



(t-t w )n(T)/2 
1 



(1 + A) 2 



4A 



v(T)/4 



ATT P (T){t-t w y+v(T)/2 

(A-l) : 



(1 + A) 2 



2T7(T)/4 



4A 



1 - 



2(1 + A) 



(16) 



where A = t/t w . More precisely, in the time regime a 2 <C 
t — t w <S^ t w , i.e. A ~ 1, 



C(t, t v 



1 



1 



{t-t w Y+^ T )/ 2 ' 



(17) 



X(t,t w ) ~1. 



This is the quasi-equilibrium regime. The dynamic func- 
tions decay at long time separations t — t w ^>t w (A > 1) 



3 



as power laws 



1 



c(t,t w ) - tn{T}/i , 

R(t,t w ) ~£ 1+J? ( T )/ 4 ) 

X(t,t w ) =X(X) > 1. 



(18) 



The crossover between the two regimes takes place at 
time t — t w ~t w . Hence, the correlation functions at dif- 
ferent waiting times do not superpose: the system ages. 




FIG. 1. Analytical result (0) for the correlation functions 
for ordered initial conditions and different t w . The dashed 
lines are power laws with exponents — r)/2 (quasi-equilibrium 
part) and —rj/i (ageing part). We have chosen a = 1, 
r?/2 = 0.03, and t w = 1, 3, 10, • • ■ , 3000, from top to bot- 
tom. 




0.85 0.9 

C(t, t u 

FIG. 2. Analytical result for the parametric Suscepti- 
bility/Correlation plot, with a = 1, 7//2 = 0.03 and 
t w = 1, 3, 10, 30 and 100, from top to bottom. The dashed 
line is the equilibrium FDT between \ and C. 



This scaling is in complete agreement with the results 
obtained by Godreche and Luck [fL2| for the spherical 



model at the critical point. Its physical interpretation is 
straightforward. At time t Wl fluctuations of wavelength 
smaller than £,{t w ) have equilibrated, while others have 
not. The quasi-equilibrium regime corresponds then to 
equilibrium fluctuations of small wavelengths. At equi- 
librium, the correlation decays as C eq (t) ~ £-2/3/^ p or 
the 2D XY model ||], one has 0/u = r?(T)/2, z = 2; 
hence the exponent describing the stationary regime of 
the correlation function is —2(3/vz = — ?y(T)/2, as found 
in Eq. ©. 

At the crossover time t — t w ~ t w , the autocorrelation 
has the value C(t w + t Wl t w ) ~ t w ri ^ T ^ 2 ■ It is interest- 
ing to compare this result with the case of ferromagnets 
quenched below T c where the crossover takes place for 
C(t w +t w ,t w ) r~j m cq (T) 2 . In this case the correlation 
function displays a plateau at the value qEA = "^cq(T) 2 , 
where qEA is the Edwards- Anderson parameter. There 
is then a clear time separation between the fast ther- 
mal fluctuations of very short wavelength (recall that 
£eq(T) *C i 1 / 2 ) and the slow relaxation due to the do- 
main wall motion. For T = T c , there is no such distinc- 
tion but only a 'pseudo-plateau' in the correlation, in the 
sense that the crossover takes place at a waiting time de- 
pendent value of the correlation. This reflects the fact 
that the order parameter vanishes in the thermodynamic 
limit. 

The last regime, which takes place at a i^-dependent 
time scale is the ageing regime. It is characterized by the 
loss of time-translation invariance and a breakdown of the 
FDT, X(t, t w ) ^ 1. On this time scale, the decorrelation 
results from the growth of the correlation length. Note 
finally that although the response function is the product 
of a stationary and an ageing part, the latter is a very 
slowly varying function of the variable t/t w . In a plot of 
X(t, t w ) vs t — t w , there would be apparently no ageing 
in the response function. 

For comparison with the numerical results of the full 
XY model, we show in Fig. |l| the analytical result ( fTl| ) 
for the correlation functions, while in Fig. H we build 
the usual plot to study the FDT violations a. We have 
plotted the susceptibility x(t, t w ) = J t dt'R(t,t') as a 
function of C(t, t w ) parameterized by the time difference 
t — t w , for different t w . With the generalized form of the 
FDT, Eq. (0), the relation between \ and C becomes 



1 r* 

x(t,t w ) =- j dt'X(t,t') 



or 



T 



(19) 



dC'x(C), 



C(t,t„) 



where the second line holds when the FDR becomes a 
function of C only. In this last case, the parametric plot 
reaches a master curve x(C) independent of t w , whose 
derivative d\/dC is given by x{C). At equilibrium, x = 1 
and x(C) = (1 - C)/T. When the FDR is not a single 
variable function, only the first line of Eq. ( pj| ) can be 
used and the FDR X(t, t w ) is not the slope of the para- 



4 



metric plot. It is clear from Eqs. ( |16[ ) that this happens 
in the 2D XY model, and a master curve x is n °t reached. 

Finally, we comment on the rather peculiar fact that 
in the time regime A 3> 1 the FDR is larger than 1, which 
means from Eq. ( Jigj ) that the susceptibility is larger than 
its equilibrium value, as can be seen in Fig. |^. When the 
system is quenched, all degrees of freedom which were 
equilibrated at temperature T = 0, fall out of equilibrium 
at time t = 0. As t increases, 'fast' degrees of freedom 
equilibrate, while the 'slow' ones are still out of equilib- 
rium. According to the interpretation of the FDR as an 
effective temperature |36| , x > 1 means that the slow de- 
grees of freedom are 'quasi-equilibrated' at a temperature 
T e ff = T/x < T. Applied to a quench from a high tem- 
perature state, this reasoning implies a FDR smaller than 
one, and the recovery of the usual case x < 1. This argu- 
ment is only qualitative, since the interpretation of FDT 
violations in terms of an effective temperature relies on 
the existence of well separated time scales, each one asso- 
ciated with an effective temperature; here instead there 
is no clear distinction between 'slow' and 'fast' degrees 
of freedom, since there are fluctuations of every length 
scales, each having its own relaxation time. 



III. AGEING DYNAMICS: MONTE-CARLO 
RESULTS 

In this section we present the results of Monte-Carlo 
simulations for the full XY model which allows to go 
beyond the spin-wave approximation. We first consider 
a quench from a completely ordered initial condition to 
check the validity of the scalings obtained in the previous 
Section. Then we consider a quench from a random initial 
condition in order to see how the vortex contribution 
changes the dynamical scaling and the FDT violations. 



A. Model and details of the simulation 

The two-dimensional XY model is defined by the 
Hamiltonian 



(20) 



where the sum runs over nearest neighbors of a square lat- 
tice of linear size L, and the spins {0 ■, j = 1, 2, . . . , L 2 } 
are bidimensional vectors of unit length. Introducing an- 
gular variables through ■ = exp(i6j), the Hamiltonian 
becomes 



H 



(i,3) 



(21) 



The Kosterlitz-Thouless transition temperature is lo- 
cated at T KT ^ 0.89 0. The model (|l|) is simulated by 
the following Monte Carlo algorithm: the angular vari- 
able 0i associated to a randomly chosen site i is updated 



to a new value 9[ randomly chosen in the interval [—w, n] , 
with probability p = min {l,e-P AB }, where AE is the 
energy variation between the two configurations. We 
present results for a system of linear size L = 512. We 
carefully checked that we are in a regime with <C L. 

We shall be interested in the spin-spin autocorrela- 
tion function C{t,t w ) and in its conjugated integrated 
response function \{t, t w ). The former is easily computed 
from 



i=l 



(22) 



where the brackets denote an average over the realiza- 
tions of the thermal noise. The computation of the lat- 
ter is standard ||. At time t w , a perturbation AH = 
— JY hi ■ (fii is added to the Hamiltonian. The magnetic 
field is random: its two components are independently 
drawn from a bimodal distribution ±h. In the linear re- 
sponse regime, the integrated response function reads 



X(t,t w ) 



1 



h 2 L 2 



(23) 



i=l 



where the overline means an average over the realizations 
of the random magnetic field. We checked that the value 
h = 0.04 was small enough for the response to be linear. 



B. Ordered initial conditions 

We investigate first the dynamic behaviour of the 
model (21) starting from a completely ordered initial 
condition. At t = 0, the temperature is switched to 
T < Tkt- The evolution of the system may be followed 
in a series of snapshots in Fig. ||. At short times, the 
system is very ordered, and there is a single color in the 
first snapshot. As t increases, fluctuations appear, and 
there is clearly a growing length scale in the system. 




FIG. 3. Snapshot of the domain growth in the XY model 



at T = 0.3 for times t n 



lCT, lCT, 10 4 and ordered initial 



conditions. The system size is 512 x 512. 
encodes the variable cos(9(x,t)). 



The grey scale 
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FIG. 4. Autocorrelation function for ordered initial condi- 
tions, for waiting times t w = 10, 30, 100, 3000, 10000 (from 
top to bottom). The final temperature is T = 0.3. The 
full lines are power laws with exponent —T)/2 — —0.03 and 
—17/4 = -0.015. 
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FIG. 5. Rescaling of the autocorrelation functions of Fig. ^[ 
The dashed lines are power laws with exponents and 0.015. 
The rescaling is good, except for the smallest waiting time. 

In Fig. |J, the autocorrelation C(t, t w ) is shown for the 
temperature T = 0.3 and different waiting times t w . The 
overall behaviour of C(t, t w ) is in nice agreement with the 
analytical results obtained in the spin-wave approxima- 
tion, recall Fig. |l|. As described in the previous section, 
there are two different regimes, depending on the time 
scale. For small time separation t — t w <C t w , C(t,t w ) 
exhibits a stationary part, which is well represented by a 
power law 



C(t, t w ) 



(24) 



while at long time, the different curves do not superpose: 
the system ages. We test the scaling ( |lT[ ) in Fig. ||, where 
we plot (t — t w ) 03 C(t,t w ) as a function of the scaling 
variable t/t w . We have actually obtained the exponent 



0.03 as the value insuring the best collapse of these data. 
This value of the exponent can be compared to the self- 
consistent computation of Ref. |35| which gives i](T — 
0.3) ~ 0.026. It is in good agreement with its numerical 
estimation @. At T = 0.7, we find r)/2 = 0.085, in 
agreement with the numerical estimation of Ref. |]37f ; the 
self-consistent computation gives 0.069. 

The parametric plot obtained numerically at T = 0.3 
for ordered initial conditions is represented in Fig. ^|. It 
is also very similar to the analytical results of Fig. ||, 
confirming in particular the unusual feature of a FDR 
larger than 1. 

The conclusion of the numerical study of ordered initial 
conditions is that the off-equilibrium dynamics is satis- 
factorily described by the spin- wave approximation. This 
is no longer true with random initial configurations on 
which we focus now. 
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FIG. 6. Parametric Susceptibility/Correlation plot for or- 
dered initial conditions. The symbols are the same as in 
Fig. I 



C. Random initial conditions 

The non-equilibrium dynamics for random initial con- 
ditions turns out to be quite different. However, our pre- 
vious analytical results can be used to understand the 
scaling behaviour of the dynamic functions in that case 
too. 




« 

FIG. 7. Snapshot of the domain growth in the XY model 
at T — 0.3 for times t w = 10, 10 2 , 10 3 and random initial 
conditions. The system size is 512 x 512. The grey scale 
encodes the variable cos(9(x,t)). 
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In Fig. we present a series of snapshots for different 
times. A simple look at Figs. and || makes the differ- 
ence between the two situations very clear. The snap- 
shots of Fig. ^ are moreover different from the 'Schlieren 
pattern s' pre sented in T = simulations of the 2D XY 
model [ [27|]2S|]3l|l and resulting from the presence of vor- 
tices. In the present work, vortices are also present, but 
they are blurred by thermal fluctuations, and do not 
clearly appear in the snapshots, except at T close to zero. 

In Fig. |, we present the correlation function C(t,t w ) 
at T = 0.3 and different waiting times t w . The two dis- 
tinct behaviours, namely a stationary part and an ageing 
part, are still present. The stationary part, correspond- 
ing to the equilibrium fluctuations, remains unchanged 
and can again be represented by the power law (0) . 




FIG. 8. Correlation functions for random initial conditions 
at T = 0.3 and waiting times t w = 10, 30, 100, 300, 1000, 
3000 and 10000 from bottom to top. 
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FIG. 9. Rescaled correlation functions for random initial 
conditions at T = 0.3; the data are the same as in Fig. ^| 
The two dashed lines are power laws with exponents and 
-0.54. 



The ageing part is different because of the presence 



of vortices in the system. By analogy with the previous 
case, we conjecture that C{t,t w ) has the scaling form 



C(i, t w ) 



1 



at) 



<I> 



(25) 



where $ is a scaling function and is the growing cor- 
relation length. Previous numerical and analytical stud- 
ies have shown that an effect of the vortices is to change 
the growth law for the correlation length from the power 
lawf(t) ~ t 1/z to the form^(i) - y/t/hxt ||M|l|. The 
scaling prediction for C(t, t w ) is tested in Fig.lTwhere we 
plot (t— t w )°- 03 C(t, t w ) as a function of includ- 
ing the logarithmic correction due to the vortices. The 
rescaling is very good. It is seen in particular that the 
long time decay is well described by a power law, i.e. 
<J>(x) - x- x for large x. At T = 0.3, we find A = 0.54. 
An asymptotic power law behaviour was also found in 
Rcf. [jl2| for the ferromagnetic spherical model. 

It is also interesting to note that the i/i^-scaling, usu- 
ally called 'simple ageing' and encountered for an 
ordered initial state, is changed to a 'non-simple' age- 
ing with the scaling tlnt w /t\nt when topological defects 
(here vortices) are present. Although the two scalings 
are equivalent in the large waiting time limit, the lat- 
ter is preasymptotically equivalent to a sub-ageing be- 
haviour |B8fl. 




0.4 0.6 

C{t,t w ) 

FIG. 10. Parametric Susceptibility/correlation plot for 
T — 0.3 and random initial conditions. The symbols are the 
same as in Fig. ^. The two full lines are similar data obtained 
with the 3D Edwards- Anderson model (see text). 



We turn now to the FDT violations. We present in 
Fig. |l^ the parametric plot Susceptibility/Correlation. 
There are still two distinct behaviours depending on the 
time scale. At short time separation FDT holds, whereas 
at long time separation a non-equilibrium regime is en- 
tered. This plot is very different from the one obtained 
in Fig. ^ for ordered initial conditions in the sense that 
the susceptibility is smaller than the value obtained from 
the FDT. In short, this means that X(t,t w ) < 1. In the 
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qualitative picture developed above in terms of an ef- 
fective temperature, this means that the fluctuations of 
wavelength larger than £(t w ) are 'quasi-equilibrated' at a 
temperature higher than the thermal bath temperature, 
keeping then the memory of their T = oo initial state. 

It is interesting to note the close similarity of this plot 
with those observed in simulations of the 3D Edwards- 
Anderson model ||,||^,[hi. We have superposed as con- 
tinuous lines the data obtained in Ref. 1(4(1 . 



IV. THE LINK BETWEEN STATICS AND 
DYNAMICS 



scaling form 



Q(m,L) = a- x T 



(30) 



where Tlx) is a scaling function whose functional from 
is well approximated by a generalized Gumbel distribu- 
tion 



Tlx) = K{e y - eV Y /2 , y = b(x-s), 



(31) 



where K, b and s are numerical constants ensuring the 
normalization of the distribution. 



In this Section we compute the finite-size Parisi func- 
tion and compare it with the finite-time fluctuation- 
dissipation ratio obtained in the previous Section. 



A. The Parisi function in the finite-size 2D XY 
model 

In order to compute the Parisi function P(q, L) for the 
2D XY model, one has to carefully define the overlap q. 
Two different overlaps between configurations (1) and (2) 
may be defined, namely 



i=l 



(26) 




or 



FIG. 11. The Parisi function P(q, L) for the 2D XY model 
at T = 0.3 and system sizes L = 4, 8, 16, 32 from right to 
left. 



q 12 = ^ E C0S ^ - ^ - 9 i + ^ > ( 2? ) 



where ip a = N -1 J2i According to the prescription of 
Ref. ||, the 0(2) symmetry of the Hamiltonian has to be 
taken into account. We have therefore considered the first 
definition. The Parisi function P(q, L) can be related to 
the probability distribution function of the magnetization 
Q(m, L), where 



1 N 

-Y 



cos 



% - i>] ■ 



(28) 



The function Q(m, L) gives the magnetization fluctua- 
tions in the direction defined by the mean orientation ip, 
which contains the physics of the critical fluctuations |^3) . 
It has been analytically computed in Refs. (23). The two 
first moments (m) and a = y/ (to 2 ) of the distribution 
scale as |B4f 



-v(T)/2 



Tim) 



(29) 



The temperature and size dependence of the distribution 
enters only these two moments, and Q(m, L) takes the 



From Q(m, L), it is easy to obtain the critical fluctua- 
tions of the overlap in the direction defined by the mean 
angle ip, i.e. using the definition by the relation 



(32) 



This function is displayed in Fig. [Tl] for different sys- 
tem sizes. Its scaling properties follow from those of 
Q(m, L) 23 §3] and are typical of the fluctuations of the 
order parameter at a critical point [fL8| . The mean val- 
ues scale with system size as (q) ~ L~ 2/3 / u , as do the 
widths of the curves. This result is expected from the 
hyperscaling relation between critical exponents and is 
compatible with a diverging susceptibility. However it is 
not sufficient to give a finite amplitude for the probabil- 
ity at q = 0. As could be estimated from Fig. [ll], the 
probability function has an exponential tail for q values 
below the mean, leading to much larger fluctuations than 
for a Gaussian function, but the tail of the distribution 
is still immeasurably small at q = 0. 

When the definition (^6|) of the overlap is used, a fi- 
nite probability is generated at q = and a distribution 
reminiscent of the Parisi function for spin glasses is ob- 
tained [E0|. In this case the 0(2) symmetry is included 
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and P(q, L) takes account of the fact that the magneti- 
zation vector diffuses around the perimeter of a circle. It 
is obtained from Q(m, L) by the relation 



2 f 1 
P(q,L) = - / d 



Q(m, L) 
s/ m 4 — 



(33) 



The result is shown in Fig. [l^ for different system sizes. 
That the tail is the product of this diffusive motion and 
is not the result of critical fluctuations can be seen by 
considering the overlap between two rigid rotors of unit 
length diffusing on a circle, for which q = | cos#|, q > 
and where 6 is the angle between the two rod directions: 
this amounts to computing P(q, L) at T = 0. In this 
case, P(q, L) is proportional to the density of states of 
cos 9 on the circle and it follows easily that 



2 1 



(34) 



The intercept at q = in Fig. |l2| is quite close to the 
above value of 2/ir ~ 0.64, while the essential singularity 
at q = 1 is removed by the fluctuations in length of the 
magnetization vector. The critical behaviour manifests 
itself again in the width and position of the maximum 
probability, which scale with system size as in Fig. [n|. 
Apart from this scaling of the peak, one would expect 
the same qualitative behaviour for the low temperature 
phase of the 3D XY or 4D XY models, where symmetry 
is broken in the thermodynamic limit. 




FIG. 12. The Parisi function P(q, L) for the 2D XY model 
at T = 0.3 and system sizes L — 4, 8, 16, 32 from right to 
left. 



While this explanation of the tail in P(q,L) seems 
rather trivial it is perhaps worth pointing out that the 
breaking of symmetry in an XY model provides a very 
simple analog for the breaking of replica symmetry in a 
mean field spin glass: the XY magnet contains an infin- 
ity of pure states a which are the points on the circle. 
On heating and cooling in zero field the system sam- 
ples different points on the circle. It is true that these 



pure states are connected by a line of constant free en- 
ergy around the circle and so in principle the system can 
change pure states without jumping free-energy barriers. 
However, symmetry breaking means that the diffusion 
constant around the circle goes to zero in the thermody- 
namic limit and the system never leaves its chosen pure 
state. Moving from one pure state to another across the 
circle does involve jumping a free-energy barrier in com- 
plete analogy with the spin glass. The overlap function 
between the pure ground states of the XY ferromagnet 
is therefore given by the expression (|34|) above. We note 
finally that distances between pure states of the 2D XY 
model on the circle do not verify the ultramctric inequal- 
ity and hence are not organized in the hierarchical 
way they are in a mean-field spin glass. 



B. Finite-size statics vs finite-time dynamics 

We are now in position to test the relation (|^) 
presented in the Introduction. We compare the dy- 
namic parametric plot Susceptibility/Correlation (data 
of Fig. [l(J) at two different temperatures, T = 0.7 and 
T = 0.3, with the curves obtained through the relation 




FIG. 13. Comparison between the dynamical data and the 
curves built from the Parisi functions using the conjecture (^j) 
at T = 0.7. The linespoints are the dynamical data, while the 
full lines are the static data. 



At T = 0.7, we find an excellent agreement between the 
dynamic and static data (see Fig. |l3| ). This strongly sup- 
ports our conjecture that the dynamic behaviour can be 
described in terms of the static distribution of a finite-size 
system of linear dimension such that L = £(t w ). There 
is a small systematic deviation between the two sets of 
curves on the left hand side of the parametric plot. This 
corresponds to the limit where the total simulation time 
is a lot greater than t w , which is most easily seen when t w 
is small. In the regime where C(t, t w ) is small and hence 







S> £(t w ), one can no longer expect the relation (g) 
to be satisfied. 

At much lower temperature T = 0.3, the two sets of 
curve coincide in the quasi-cquilibrium regime C > c/ea 
and the departure point from the FDT is correctly pre- 
dicted by the static data (see Fig. |l4|). As discussed 
above, this breaking point is given by C(t w +t w ,t w ) ~ 

t w , which coincides with the peak in the Parisi func- 
tion (q) = L~ 2f} / U , with L = tlJ z . Although the qualita- 
tive shape of the curves in the ageing regime is the same 
for the two sets, there is a quantitative discrepancy. This 
suggest that, while Eq. (||) is a correct starting point, the 
analytic P(q, L) that we have calculated is not in this case 
quantitatively accurate. 




0.4 0.6 

C(t, t u 

FIG. 14. Comparison between the dynamical data and the 
curves built from the Parisi functions using the conjecture ([]) 
at T = 0.3. The linespoints are the dynamical data, while the 
full lines are the static data. 



Indeed, in deriving this function we have explicitly ne- 
glected vortices which play a crucial role in the transition 
from order to disorder in the 2D XY model |l9|. Since 
vortices are present during the coarsening process, one 
should consider their contribution to the static fluctua- 
tions. We do not know precisely how this has to be done, 
for instance we do not know what density and what vor- 
tex correlation to introduce. However, it is clear that the 
presence of defects can dramatically alter the shape of the 
function P(q,L). As an illustration, we have calculated 
numerically P(q, L) in the presence of a pinned defect. 
We have slaved a single spin in the opposite direction to 
its left neighbour and found that the tail at small q is 
considerably reduced. 

As a first step in investigating the effect of vortices, we 
show in Fig. 0| vortex configurations at the same waiting 
time, t w = 1(F\ for T = 0.3 and T = 0.7. While the num- 
ber of vortices is of the same order (474 for T = 0.3 and 
602 for T — 0.7), their correlations are clearly different. 
At low temperature, the vortices of opposite charge are 
essentially uncorrelated, while at higher temperature, the 
majority of vortices are in closely bound pairs. This is 



because the diffusion of vortices is more efficient at higher 
temperature and so their natural tendency to form pairs 
is largely satisfied. The effect of a free vortex in the 
system is much more dramatic than a vortex pair [^9| . 
In Fig. |l5|, there are far more free vortices at the lower 
temperature and hence one can expect that the static 
analysis in the absence of vortices describes more accu- 
rately the situation at T = 0.7 than at T = 0.3. This is 
consistent with our observations in Figs. [l3] and |l4| . 
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FIG. 15. Vortex configurations at time t w = 10 after 
a quench from a random initial condition to a temperature 
T — 0.3 (above) and T = 0.7 (below). Antivortices are repre- 
sented by squares, vortices by plus symbols. 



V. DISCUSSION 

We have studied in this work the nonequilibrium crit- 
ical dynamics of the 2D XY model and its relation with 
static properties. We have investigated the ageing be- 
haviour of the system by looking at two-time dynam- 
ical functions and by considering quenches from ran- 
dom and ordered initial states to temperatures below the 
Kosterlitz-Thouless transition. We have then shown that 
the finite-time violations of the FDT are related to the 
finite-size equilibrium Parisi function. We discuss now 
some consequences of our results. 

A first important feature is the scaling (^5|) of the spin- 
spin autocorrelation function. Figure @ shows that the 
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behaviour of C(t,t w ) is, at first sight, very similar to 
the usual additive decomposition between a stationary 
and an ageing part M, We have seen that, at a crit- 
ical point, this decomposition is instead multiplicative. 
Interestingly, the same scaling is observed in numerical 
simulations of the 3D Edwards- Anderson model [§J],fl3| . 
It is indeed notoriously difficult to obtain the value of 
the Edwards- Anderson parameter qEA, the plateau in the 
correlation function, in 3D Ising spin glasses p2| . This 
similarity means that the dynamics of 3D spin glasses is 
influenced, at least at short time, by the presence of a 
critical point. 

This multiplicative decomposition has also conse- 
quences in the scaling behaviour of the response func- 
tion. In the case of a quench in a ferromagnetic phase, 
the susceptibility x(A t w ) can be decomposed into a 'bulk' 
and an 'interface' contribution [||. The bulk contribution 
accounts for the equilibrium fluctuations inside the do- 
mains, while the interface contribution accounts for the 
domain wall response. Since their density vanishes at 
long time, so does their response function. This is not 
true at the critical point, where there are no clearly well 
defined 'domains': on a length scale the 'bulk' is 
composed of critical fluctuations, a self-similar structure 
composed of domains within domains Q|l5j[l^] . The 
direct consequence is that in a parametric Susceptibil- 
ity/Correlation plot, no master curve x(C) ( as predicted 
by the dynamical mean-field theory) can ever be reached. 
The difference is irrelevant for the Ising spin chain at 
T = 0, since there is no stationary part at all fl^ , |l3|| . 

We have conjectured a finite-time, finite-size general- 
ization of the relation (||) in Eq. (||) and shown that it 
works in the present model provided that the tempera- 
ture is not very low. We wish to mention that the con- 
jecture (J3j) has been implicitly used by Marinari et al. in 
Rcf . ||] . These authors have indeed shown that dynamic 
and static data for a 3D spin glass perfectly overlap, us- 
ing the relation (||) with t w = 10 5 and L = 16. Since 
P(q, L) still depends on the size L for L = 16, and the \ 
vs C plot still depends on t w at t w = 10 5 (because of the 
scaling (p5|)), the coincidence of the two sets of curves 
necessarily results from Eq. (||). It is interesting to note 
that this correspondance provides a direct measure of the 
dynamic correlation length 

As the temperature decreases, the agreement between 
the static and dynamic data is not as good because the 
density of free vortices increases. This is consistent with 
the results for the Ising chain at T = jl3| , where Eq. (||) 
has been used to deduce from the dynamic FDT viola- 
tions a static 'Parisi function'. The latter is not linked to 
the T = fluctuations of the spin glass order parameter 
which are trivial f44f| . Rather, the dynamic behaviour 
has been shown to be entirely controlled by the domain 
walls. It would be interesting to see if a similar feature 
occurs at low temperature in the 3D Edwards- Anderson 
model. 

Since the ageing dynamics of the 2D XY model closely 
resembles that of more complex systems, it is important 



to understand why this is so. We believe that the main 
ingredient for this behaviour is the presence of the crit- 
ical fluctuations. After a quench, magnetic fluctuations 
of wave vectors between k ~ 1/a and k ~ 1/L try to 
develop in a system of size L and lattice spacing a. At 
time t w , only the fluctuations with k > t w ^^ z are equi- 
librated, while those with k < t w are out of equilib- 
rium, the smallest k being very far away from the equilib- 
rium. As a result, the decay of the autocorrelation func- 
tion is dictated by the growth of the correlation length 
which gives the )-scaling of Eq. (p5|). Hence, 

this decay takes place in a single time scale. However, 
the coexistence of multiple length scales, each having its 
own relaxation time, implies FDT violations described 
by the smooth parametric plots of Fig. The coexis- 
tence of a single time scale in the ageing regime together 
with a smooth \ vs C pl°t arises naturally within this 
critical scenario. This feature, observed in 3D and 4D 

gum 



spin glasses 



has no interpretation within 
the mean-field dynamical theory, where multiple effec- 
tive temperatures are necessarily associated with multi- 
ple time scales @,|,|,|o|. 

The relevance of multiple length scales in the age- 
ing dynamics of spin glasses has been emphasized by 
Bouchaud Ej|. The example of the 2D XY model has 
shown that these length scales account well for the scal- 
ing properties of two-time dynamical functions and FDT 
violations. However, a necessary ingredient to explain 
more refined experimental results, such as the 'rejuve- 
nation and memory' effects |^(J, lies in the role played 
by the temperature [45] ]. In this picture, based on the 
behaviour of elastic objects pinned by disorder, the re- 
laxation time t{£) associated with a length scale I scales 
as t(£) cx exp [T(T)£ e /T], in the notations of Ref. ||. 
In the case of the 2D XY model the scaling is instead 
t(t) ~ £ z , independently of the temperature T. This 
means that such a model will not be able to capture sub- 
tle temperature-cycling effects. In this respect, it would 
be interesting to combine the modified domain growth 
approach recently proposed by Yoshino et al. [|47| with 
nonequilibrium critical dynamics. 
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